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Abstract 


For over 25 years, NASA has supported the development of space 
radiation transport models for shielding applications. The NASA space 
radiation transport model now predicts dose and dose equivalent in 
Earth and Mars orbit to an accuracy of ±20%. However, because larger 
errors may occur in particle fluence predictions, there is interest in 
further assessments and improvements in NASA’s space radiation 
transport model. In this paper, we consider the effects of the isotopic 
composition of the primary galactic cosmic rays (GCR) and the isotopic 
dependence of nuclear fragmentation cross-sections on the solution to 
transport models used for shielding studies. Satellite measurements are 
used to describe the isotopic composition of the GCR. Using NASA’s 
quantum multiple-scattering theory of nuclear fragmentation 
(QMSFRG) and high-charge and energy (HZETRN) transport code, we 
study the effect of the isotopic dependence of the primary GCR compo- 
sition and secondary nuclei on shielding calculations. The QMSFRG is 
shown to accurately describe the iso-spin dependence of nuclear frag- 
mentation. The principal finding of this study is that large errors 
(+100%) will occur in the mass-fluence spectra when comparing trans- 
port models that use a complete isotope grid (~1 70 ions) to ones that 
use a reduced isotope grid, for example the 59 ion-grid used in the 
HZETRN code in the past, however less significant errors (<20%) occur 
in the elemental-fluence spectra. Because a complete isotope grid is 
readily handled on small computer workstations and is needed for 
several applications studying GCR propagation and scattering, it is 
recommended that they be used for future GCR studies. 


Introduction 


An 


important goal for NASA’s Space Radiation Health Program 
is to develop a predictive capability to predict the galactic cosmic 
rays (GCR) fluence spectra to within a +25% accuracy (Anon., 1998). NASA has developed the 
HZETRN (high-charge and energy transport) code (Wilson, 1977; Wilson and Badavi, 1986; 
Wilson et al., 1991, 1995a) as a science application and engineering design tool (Wilson et al., 
1993) to be used in space radiation shielding studies. HZETRN has been validated in its ability 
to predict total dose and dose equivalent behind several materials in space to within +20% on 
multiple space missions in Earth orbit (Table 1) and the MARIE experiment on the Odyssey 
spacecraft in Mars orbit (Figure 1) (Cucinotta et al., 2000a, Badhwar and Cucinotta, 2000, 
Badhwar et al., 2001, Zeitlin, et al., 2002, http://srhp.jsc.nasa.gov). However, interest in 
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fluence-based approaches to risk assessment (Cucinotta et al., 1995a, 1996a) suggests that more 
stringent tests of transport code accuracy need to be made, and the quantities dose and dose 
equivalent are deemed as necessary, but not sufficient tests of their accuracy. In this regard, we 
note that dose and dose equivalent are integral quantities that receive contributions from many 
GCR charge groups. Large uncertainties currently exist in radiation quality factors (Cucinotta, 
et al., 2001) and current methodologies to estimate health risks such that dose and dose equiva- 
lent may be insufficient as tests of transport code accuracy. The use of ion fluence as a basis for 
tests for accuracy provides for sufficient generality to ensure accuracy in GCR transport models, 
including under the circumstances of revision of radiation quality factors or integration of 
alternative risk assessment approaches in the future. 


MARE Measurements and Model Calculations 



Figure 1: Comparison of HZETRN code predictions to MARIE experiment (Zeitlin eta/., 2003) on 
Odyssey orbiting Mars. 

In the description of the transport of the GCR in shielding materials or tissue, a common ap- 
proximation used in the past is to consider only the elemental composition of the primary GCR 
and a reduced isotope grid for the secondary nuclei produced in nuclear fragmentation. In this 
paper, we analyze the role of the isotopic dependence of the GCR primary composition and 
nuclear fragmentation in predicting the fluence of the GCR behind arbitrary shielding configura- 
tions. Our study is an important milestone in achieving NASA’s goal of accurate GCR transport 
codes since, for the first time, a complete isotope grid has been achieved in a GCR transport 
model and we document the error inherent in former approaches. Also, for applications that will 
consider radioactive isotopes produced in the atmosphere or shielding, our study provides a 
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useful tool to perform such analyses. Other applications where non-stable nuclei are considered 
are studies of the origin of the GCR where so-called cosmic-ray “clocks” consider the primary 
or secondary GCR with life-times on the order of confinement time in the galaxy (~1 M-yrs) 
(Yanasak et al., 1999). Several GCR “clock” nuclei, including l0 Be and 26 Al, were not consid- 
ered in the isotope grid used in HZETRN in the past. Finally, the MARIE experiment orbiting 
Mars on the Odyssey spacecraft (Badhwar, 2003, and Zeitlin et al., 2003) is collecting new data 
on the GCR environment near Mars, and our study provides an opportunity to begin new inves- 
tigations on the accuracy of computational models used to describe the GCR environment. 

Historically, the HZETRN code grew from a 29-ion isotope grid used in the 1980s and early 
1990s (Wilson, et al., 1991) to an extension to a 32-ion isotope grid made in 1993 in order to 
include all light ions (Cucinotta, 1993). Because of the limitations of random access memory 
(RAM) present on the computer workstations of the early 1990s, sensitivity studies were made 
for mono-energetic ion beams to study the minimum number of isotopes for convergence, 
resulting in the use of 59-isotope grid (Kim, et al., 1994), and all GCR studies since 1994 have 
used the 59-isotope grid (Shinn, et al., 1994). However, there are several reasons to reconsider 
the use of the full isotope grid for GCR transport problems. First, the isotopic dependence of 
the primary GCR has not been considered in past shielding studies and may lead to errors in the 
description of both primary ion attenuation and secondary particle production, including the role 
of high-energy neutron production from the many neutron-rich species that occur. Secondly, the 
studies of Kim et al. (1994) used the NUCFRG2 model of fragmentation (Wilson et al., 1994), 
which does not provide a correct description of the even-odd effect observed in fragment pro- 
duction or of the projectile iso-spin dependence observed experimentally (Knott et al., 1996, 
1997; Zeitlin et al., 2001). Thirdly, the sensitivity studies made by Kim et al. (1994) used a 
“calibration” of the isotope grid to 56 Fe beams; however a larger isotope grid occurs when all 
GCR projectile nuclei are considered. Fourth, the error in the range-energy and stopping powers 
that results from the use of a reduced isotope grid, although expected to be small for large mass 
numbers, A»l, is an unnecessary one for transport calculations. Finally, the improved compu- 
tational speed and memory available on current small computer workstations makes including a 
complete isotope grid in the HZETRN code possible. 

In this paper, we describe the implementation of the HZETRN code to include the full isotopic 
dependence of the primary GCR. We describe the physics of isotopic effects in GCR transport, 
and the fragmentation parameters are a key component of this description. The quantum mul- 
tiple scattering theory of nuclear fragmentation (QMSFRG) is used as the generator for frag- 
mentation cross-sections used in our study. An empirical model of the isotopic composition of 
the primary GCR, including its solar modulation, is also described. For GCR problems, we 
identified an isotope grid of 170 ions and made comparisons to previous HZETRN results using 
the reduced-grid of 59 isotopes. We note that a preliminary version of this report used an iso- 
tope grid of 141 ions (Cucinotta, 2002), which has been expanded in this current work. The 
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present code includes all of the abundant nuclei in the GCR environment with fluxes greater 
than about 1 0 2 /cm 2 /yr and nuclei produced in fragmentation events with production cross- 
sections greater than about 1 mb. Several nuclei with smaller primary abundances or production 
cross-sections, which are of interest for scientific reasons are also included in the expanded 
HZETRN model described herein. The resultant code includes many neutron-rich nuclei that 
have been ignored in the past with iso-spin components ranging from T= +3/2 to T=-3. 

Nasa currently uses the GCR model of Badhwar and O’Neill 
(1992) to describe the elemental composition and energy spectra of 
the GCR, including their modulation by the Sun’s magnetic field. 

In this GCR representation, we considered only the most abundant 
GCR nuclei for each element and counted other isotopes of identi- 
cal charge as the abundant isotope. However, theoretical models and satellite measurements of 
the GCR have long considered the isotopic composition of the GCR and their modification 
through transport in interstellar space, including estimating the primary nuclear composition at 
stellar sources (Parker, 1965, Webber et al., 1990a, Fields et al., 1994). The approach used here 
is to estimate an energy-independent isotopic fraction, f from satellite measurements, which are 
constrained to obey the sum rule 

(1) *(Z,£) = ^ //(4,Z,£) 

where the left-hand side of eq. (1) is the elemental spectra from the Badhwar and O’Neill model 
and 'Yjjfj = 1 . Equation (1) is used herein as an initial estimate of the influence of the primary 
isotopic composition on GCR shielding calculations. Experimental studies have included 
measurements on the Pioneer, Voyager, and Ulysses spacecraft. A survey of such data (Hesse et 
al., 1991; Lukasiak, et al., 1993, 1995; Webber et al., 1985, 1990; Wiedenback et al, 1981, 

1985) was made, with the results shown in Table 2. For this compilation, we note that, since 
secondary fragment production is modulated by the transit time in the heliosphere, the isotopic 
fraction is dependent on the position in the solar cycle. In Table 2 we have used data on isoto- 
pic fractions near solar maximum where the isotopic fraction for nuclei produced within the 
heliosphere are expected to be at a maximum (Hesse et al., 1991; Lukasiak et al., 1993; Webber 
et al, 1985, 1990; Wiedenback et al, 1981, 1985). For several elements listed in Table 2, 
information on solar modulation was not available and isotopic fractions for near-Earth and 
GCR sources were set identical, using the near-Earth estimate. 

We next discuss some areas of interest for understanding the GCR isotopic composition. The 
GCR path-length distribution represents the mean amount of interstellar and interplanetary 
material intersected by cosmic rays prior to their arrival in the near-Earth environment. This 
distribution is currently estimated to vary between 3 and 20 g/cm 2 of approximately 90% H and 
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1 0% He. The variability in the estimates is because they are model- dependent on nuclear frag- 
mentation parameters, cosmological interaction terms, and factors related to the near- Earth or 
deep space measurements, including energy and mass resolution of the detectors as well as 
position in the solar system. In most current models, the GCR path- length distribution is esti- 
mated to be velocity- and rigidity-dependent (i.e., energy- dependent). A common form for the 
GCR path- length distribution is (Webber et al . , 1990a) 

X = 10.8)3(— )-°- 6 for R>4GV; X =10.8/3 for R<4GV 

R o 

where /3 is velocity relative to the speed of light and R the rigidity. 

Of note is that, based on the solar modulation theory of Parker (1965), the inter-planetary por- 
tion of this distribution would be modulated over the solar cycle, suggesting the form 

A(0, /*) Xfaigrgignu,. + Xj nter pi anetar yGt), r) 

where <J> is the solar modulation parameter and r (in A.U.) is the radial distance from the Sun. 

It then follows that the isotopic abundances are dependent on the modulation parameter and 
radial distance, and will be energy- dependent. The so-called Leaky Box Model is the approach 
most often used to estimate the path- length distribution. Here the interstellar propagation of the 
GCR, including source terms, ionization and nuclear scattering terms as well as other cosmo- 
logical interaction terms, is described. However, the separation described above between 
scattering inside and outside of the heliosphere is not made. Because the path length is energy- 
dependent, the Leaky Box Model divides interstellar materials into slabs and uses weighted 
average over distinct energy bins to reconstruct the GCR abundances near Earth (or the inverse 
problem to reconstruct the source abundances). Such an approach is not easily implemented by 
the HZETRN code, but could be treated by the non-perturbative Green’s function approach to 
HZE transport of Wilson et al. (1994). Herein, these considerations are not addressed and the 
solar modulation of the GCR abundances are treated in an energy- independent manner using 
empirical estimates based on satellite data. To go beyond the present or possibly alternative 
parametric approaches will require implementing the Leaky Box Model to create the boundary 
condition for the HZETRN shielding calculations. 

Based on the results of Lukasiak et al. (1993), the following empirical formula is used to 
describe the dependence of the isotopic ratio (near- Earth) on the solar cycle including a depen- 
dence on the modulation parameter O(MV) 

(2) f(A,Z) = f nm (A,Z)+\ex f (yM-\\ 
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where y is found by interpolating from the source and near solar maximum values. The second 
term in eq. (2) can be positive or negative as dependent on the conditions for a net gain or 
reduction in the cosmic nuclei during transit from the galactic sources (Fields et al., 1994) to the 
inner heliosphere. For the Z=1 and Z=2 ions, we use the following empirical formula to esti- 
mate the primary (near-Earth) 2 H and 3 He spectra 

(за) &„,(£) = {O.0764+ 0.097exp[-0.5(ln(£/1660)/1.306) ! ]}^ 4 „,(£) 

( зб ) &„(£) = 0 . 2 &„,(£) 

where <p }He (E) and <p 2H (E) are subtracted from the Badhwar and O’Neill model (1992) for 4> 4He (E) 
and (f> IH (E), respectively. Figure 2 shows the assumed dependence of the isotope abundance 
scaling as a function on the solar modulation parameter described by eq. (3) for )3 C and 15 N. 
Similar results are found for other isotopes listed in Table 2. Examples of the GCR energy 
spectra for hydrogen and helium are shown in Figure 3a and for the isotopes of Ne, Si, and Fe 
at solar minimum in Figure 3b. 



Figure 2: Parametric model for describing the change in isotopic composition with the solar 
modulation parameter, g>(MV). 
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Figure 3a: Energy spectra for hydrogen and helium isotopes near solar minimum (<{>=428 MV) 
and solar maximum (g>=1050 MV). 
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Figure 3b: Energy spectra for Ne, Si, and Fe isotopes near solar minimum (<j>=428), showing 
contributions from different isotopes to primary GCR composition. 






T 

A. he reduction of the full 3-dimensional (3D) Boltzmann 
transport equation within the straight-ahead and continuous 
slowing down approximations and assuming velocity con- 
servative fragmentation events accurately describes the 
transport of GCR heavy ions (Wilson et al. 1986, 1991, 1995). In this model, the heavy ion 
flux, <j)/E,x j of an ion j with mass number A., charge number Z. energy E (in units of MeV/u) at 

shielding depth x (in units of g/cm 2 ) is determined by the partial differential equation (Wilson et 
al., 1991) 

where S.(E) is the stopping power of ion j, a fE) is the energy-dependent absorption cross- 
section (cm 1 ) and a ..(E) is the fragmentation cross-section for producing an ion j from k. The 
solution to eq. (4) in the HZETRN code is found using the methods of characteristics where the 
coordinate transformation 

(5) rj. = x-RjiE ); 4j = x + Rj(E) 

and the scaled flux 

(6) = SjiEWji^E) 

are introduced, leading to the transport equation 


Isotopic Effects 
in GCR Transport 


( 7 ) 



IxOij’tj) 


k v t 


where v. = Zj/A.. The numerical solution to equation (7) is implemented using a marching 
procedure as described by Wilson et al. (1990). The solution of the GCR transport equation for 
light ions and neutrons is distinct from eq. (7) because of the broad redistribution of energy in 
collision events and is described elsewhere (Wilson et al., 1991, Cucinotta et al., 1994a, and 
Clowdsley et al., 2000). 

Because eq. (7) is a coupled integro-differential equation for the many GCR primary and sec- 
ondary nuclei, required computer memory allocations increase rapidly as the number of ion 
species is increased and was an important consideration in the early 1990s. However, such 
practical limitations no longer exist, even on small computer workstations with sufficient RAM 
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(e.g. Pentium III or higher), and an unrestricted isotope grid can presently be implemented with 
no memory or storage problems. When using a reduced isotope grid, fragments not contained in 
the grid are assigned to a nearby mass of identical charge number. This introduces two types of 
errors: First, the range-energy and stopping powers are altered from their true values. This error 
is expected to be small at high-energies for A» 1 since here ion ranges are well described by 
A/Z 2 scaling factors. The second type of error occurs in the absorption and fragmentation 
cross-sections. Here, because of the reassignment of the mass number when using a reduced 
grid, an error is introduced by the change in neutron number from its true values. Such errors 
arise from modifying the atomic transport properties and nuclear fragmentation cross-sections, 
and spectra of high-energy neutrons produced in fragmentation events from their true values by 
forcing the physics onto a reduced isotope grid. This latter error is expected to be small for 
heavy target nuclei (.4 >16) because neutron production is dominated by light-particle (n, p, d, t, 
h, and a) interactions on target nuclei, but may be non-negligible for light target atoms (A<16) 
(Cucinotta et al. 1998a). Since materials with high-hydrogen content are known to be the 
optimal shielding materials, the changes in the neutron fluence due to the use of a full isotope 

We next discuss the quantum multiple scattering descrip- 
tion of heavy ion fragmentation (QMSFRG), which has been 
quite successful in describing the physics of the abrasion - 
ablation model of fragmentation and experimental data 
(Cucinotta, et al. 1994b, 1997, 1998b). The scattering ampli- 
tude for the heavy ion collision is related to the cross-section 
by the phase space of each particle that appears in the final state. In the QMSFRG theory for 
inclusive reactions in which a single fragment originating in the projectile is measured, closure 
is performed on the final target state with a momentum vector denoted used to represent 
these states. The cross-section is then given by 

(8) da = ^ Z d$ x dp F . X FI - E /)s(p, - p/ )|T fl | 

where is the relative projectile-target velocity, F* represents the pre-fragments formed in the 
projectile-target interaction, n is the number of nucleons knocked out of the projectile in the 
overlap region with the target (Cucinotta, 1994c), and i and /label the initial and final states, 
respectively. The pre-fragments decay through particle emission if sufficient energy is avail- 
able. To include the phase space of decay products of F*, we write 

(9) dp F . = dp F \[dp r 

r=0 


grid should be considered. 

Isotopic Effects 
in Quantum 
Fragmentation 
Model 
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where r is the ions (if any) emitted in the decay of the F*. In considering nucleon production 
from the decay, we would study the p r (Cucinotta et al., 1997b). The total momentum transfer 
is q = p r -p r , where p T is the initial target momentum. We use the momentum conserving 
delta-function in eq. (8) to eliminate or one of the p >, from eq. (8). The 3D inclusive cross- 
section for producing an ion F is given without approximation by 


( 10 ) 


dc_ 

dp F 


(2*) 4 

P 


x - E f»<p, -p/ w 2 


For elastic scattering or the excitation of discrete states, the relation between the transition 
matrix T {1 and the inclusive cross-sections is trivial. For fragmentation reactions, where several 
to many particles are present in the final state, the integrals in eq. (10) become intractable and 
approximations must be introduced. One approach is to use a closure approximation on all 
unobserved projectile fragments, however at the expense of losing information on final state 
interactions among the projectile fragments. Real progress in reducing the multiparticle mo- 
mentum integrals to a computationally feasible form is achieved only after studying the struc- 
ture of the nucleus-nucleus transition matrix. The equations of motion for nuclear scattering are 
expressed in terms of the transition operator, which represents an infinite series for the multiple 
scattering of constituents of the projectile and target nucleon. The strong nature of the nuclear 
force requires a non-perturbative solution to the scattering problem. A relativistic theory is of 
interest for the space radiation databases, because of the high energies of the particles and the 
large number of production processes that are naturally included in a relativistic theory. A 
relativistically covariant formulation of the problem has been put forth by Maung and co-work- 
ers (1996) using meson exchange theory. The basic approach, in both relativistic and non- 
relativistic multiple scattering theories, is to re-sum the multiple scattering series, which is 
expressed in terms of the irreducible and reducible exchange diagrams in the RMST or the 
nuclear potential in the NRMST, in terms of the transition matrix for projectile and target nuclei 
constituents. This avoids having to deal directly with the highly singular behavior of the nuclear 
potential at short distances, and instead the constituent transition matrix is used, which is known 
from experiment. The integral equation approach is quite successful for studying elastic scatter- 
ing where a one-body integral equation can be found by formulating an optical potential. For 
studying knockout and fragmentation reactions, the Eikonal approximation (Cucinotta et al., 
1989) is useful in order to reduce a many-body integral equation to a manageable form. The 
importance of final state interactions between projectile fragments suggests the use of a Faddeev 
type integral equation. 

In the RMST, the infinite sum of meson exchange diagrams is written as an integral equation of 
the Bethe-Salpeter form (Maung et al., 1996). The Bethe-Salpeter equation is reduced to a 3D 
form using a covariant 3D relativistic propagator. The propagator of Maung et al. (1996) is 
most useful for performing the 3D reduction, since it treats the target and projectile constituents 
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on an equal footing, avoiding non-physical singularities that occur with other propagators. The 
transition operator derived in the RMST is written as 

(11) T = K + KGT 

where G is the Bethe-Salpeter propagator representing the two nuclei in intermediate states and 
the kernel K is the sum of all irreducible diagrams based on meson exchange theory for scatter- 
ing of the projectile and target constituents. The kernel is decomposed into various terms 
corresponding to one meson exchange between constituents, two meson exchanges between 
constituents, two meson exchanges between more than one constituent, 

(12) K = K ] +K 2 +K 2X +... 

This infinite sum of irreducible diagrams is described by Maung et al. (1996). The 3D reduc- 
tion of the RMST is found by introducing an approximate propagator g to obtain the coupled 
integral equations: 

(13) T = V + VgT 
with 

(14) V = K + K(G-g)V 

The 3D reduction is chosen to represent the best approximation to an exact propagator G. In 
application, the approximation V = AT, is often evoked. The effects of nuclear clustering are 
considered in the MST by assuming the constituent interactions are those between clusters 
rather than the choice of nucleons (Maung et al., 1996). The RMST with clusters involves 
complicated summations over irreducible diagrams among the cluster constituents. The choice 
of which cluster configuration is chosen is determined by reaction channel and nuclear structure 
considerations. The convergence of a cluster expansion series should be more rapid than the 
nucleon one when the kernel is known; however, more detailed bound state properties may be 
involved for performing such calculations. 

The NRMST is obtained from eqs. (12)— (14) by approximating the full kernel by the leading 
order term corresponding to one-meson exchange diagrams and using a non-relativistic reduc- 
tion of the 3D propagator g. The potential term is the sum of the interactions of the constituents 

Ap A? 

os) 2X 

j= 1 0=1 

and the non-relativistic propagator is given by 
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06) g m =(E-H r -H T \' 


where H p and Hj are the projectile and target internal Hamiltonians, respectively. The constitu- 
ent interactions involve the full many-body problem as described by the integral equation 


( 1 7) — Vy + Vq g NR Ty 

where V is the nucleon-nucleon potential and the propagator includes the effects of nuclear 
binding. At high energies, the relative kinetic energy of the constituents is much larger than the 
binding energy, such that the propagator can be evaluated in the impulse approximation 

(18) g,=(E-T r -T T Y 


and the constituent interactions are replaced by the free interactions which are truly of the two- 
body form. For high-energy reactions, the scattering is often confined to the forward direction 
and the Eikonal approximation is accurate and to be used to reduce the scattering problem to a 
closed form expression. There are several approaches for deriving the Eikonal form of the 
MST. Here we continue our considerations of the nucleus-nucleus propagator and introduce the 
Eikonal propagator 


(19) 


Seik 


f k-(k-k') ) 
P« ) 


8 


(k-k') 


The insertion of the Eikonal propagator into the MST allows for a summation of the series into a 
closed form expression. An alternative coupled-channels approach to the Eikonal approxima- 
tion is considered by Cucinotta et al. (1989). Calculations using the Eikonal model are considered 
next. 


The pre-fragment excitation spectrum following nucleon or alpha particle abrasion can be 
represented in terms of an impact parameter-dependent convolution of the pre-fragment excita- 
tion response for a transition of the pre-fragment core from state n to n' and the project fireball 
response (Cucinotta et al., 1998b) 


( 20 ) 


do 


de 


F* 


= <T\jd 2 qd 2 bd 2 b' e ,w *>P,,.(b,b')A.Aq,b,b,E r .)\T > 


where b (b ’) is the impact parameter, and q the momentum transfer. The abrasion response is 
defined as the interaction of the projectile fireball with the target after performing closure over 
the final fireball states 


12 


(21) A ME,) - f^L 


< R'\Q RT + (b') \RxR\ Q RT (b ) | R' > <5(£,. -E f ) 


where the Q RT represent the fireball-target profile operator, and k R the projectile fireball momen- 
tum vector. The abrasion response represents a complicated many-body operator that is solved 
by approximation. The one-particle abrasion response has been evaluated using the shell model 
response functions (Cucinotta and Dubey, 1 993). The pre-fragment excitation is described in 
terms of the transition matrix 

(22) P ny (b,b',E f ) = <F\\Q F . T \b')\F\, xF*,\Q F . T (b)\F\ > 

where matrix elements for the pre-fragment excitation are evaluated over the many-body profile 
operators, Q ^ The model uses a convolution approach to derive the mutli-knockout spectrum 
from the single-fragmentation term (Cucinotta and Dubey, 1994b). 

The de-excitation of the pre-fragments in the QMSFRG model is described in a stochastic process 
using a Master equation for nuclear de-excitation by particle emission (Cucinotta and Wilson, 

1996b). If f(E,t) is the probability of finding the nuclei b at time t with excitation energy E b and 
P k b (E) the probability that the nuclei, b will emit ion k with energy E, then the Master equation is 

(23) = %jd E r(.Z,OP;(E)-'ZI‘IEf t (K.t)P?(E) 

In eq. (23), the first term on the right corresponds to gains by decays a -> b+j and the second 
term from losses due to decays b -> c+k where the j (or k) are light-particle emissions (n, p, d, t, 
h, or a). Eq. (23) is solved by iteration up to medium excitation energies (below 150 MeV) and 
by approximation for high excitation energies (Cucinotta et al. 1996b, 1998b). Important fea- 
tures of this solution is the correct description of the nuclear level-density including nuclear- 
shell effects at low-excitation energies, and the use of measured values for the nuclear masses. 
The fragmentation cross-section is then evaluated from eqs. (21)-(23) as 

(24) O r (A„Z r ) = 21 iE r • 

a f \z f * at F , 

where f(A F .,Z F , >A F ,Z F ) is the solution to eq. (23). 

The QMSFRG theory reduces to the OPTFRG model (Townsend et al., 1986) when energy 
conservation and nuclear medium effects are ignored and closure approximations on the pre- 
fragment and fireball states are made, and to the NUCFRG2 model (Wilson et al. 1995b) when 
the optical operators are expressed as volume overlaps functions. In both the OPTFRG and 
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NUCFRG2 models, all information on the pre-fragment excitation spectrum is lost and thus 
must be introduced in an ad-hoc manner independent of the collision model used to describe 
mass removal. Quantum interference effects in nuclear fragmentation can by ignored for heavy 
projectiles, but is important in the case of small projectile mass (Cucinotta et al., 1992a) and 
also modifies nuclear absorption in quasi-inelastic scattering (Cucinotta et al., 1 992b, 1 992c). 
For light particle fragmentation, we use the models of Cucinotta and co-workers described 
previously (Cucinotta, 1993, Cucinotta et al., 1995b). 

The iso-spin dependence in fragmentation cross-sections enters in several ways. First, the energy- 
dependent two-body nucleon interaction parameters are summed over the possible projectile and 
target nucleon scattering combinations. Second, the nuclear wave-functions will have a dependence 
on shell structure and the nuclear surface. The largest contribution to iso-spin effects occur in the 
nuclear ablation process (de-excitation) where the pre-fragments formed and their level spectra are 
greatly influenced by the projectile or pre-fragment iso-spin. Important aspects of the model include 
the use of Coulomb trajectories (Cucinotta et al., 1997a) and in-medium nucleon-nucleon interaction 
(Tripathi et al., 2000). Their effects on abrasion cross-sections are illustrated in Figure 4 for 32 S 
fragmentation on aluminum targets at kinetic energy of 100 MeV/u. Clearly these aspects of the 
physics of fragmentation are important and their inclusion leads to improvements in the accuracy of 
fragmentation and absorption cross-sections databases. 



0 2 4 6 8 10 12 


b, fm 

Figure 4: Abrasion cross-sections versus impact parameter for n=l, 2, 4, and 12 nucleon removal 
in 32 S fragmentation at 0.1 GeV/u on Al targets showing the corrections for Coulomb trajectories 
and the in-medium nucleon-nucleon interaction. 
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Results 


We first illustrate the accuracy of the QMSFRG model and the 
effects of iso-spin on fragmentation cross-sections. Figures 5-16 show 
comparisons of the model to experimental data for the elemental distri- 
butions of fragments for several nuclei of similar mass number. All comparisons include the 
one- and two-nucleon removal cross-sections by electromagnetic dissociation (Norbury et al., 
1988). The iso-spin, T=0 nuclei display large odd-even effects, which are reduced for the 
T > 0 nuclei. The odd-even effects are present for all target nuclei, however are reduced for 
hydrogen targets due to the small abrasion probability for large mass removal on hydrogen. To 
see one of the errors that results from transporting ions using a reduced mass-grid, we compare 
fragmentation cross-sections for nearby projectiles where large differences in many of the 
production cross-sections occur for neighboring projectile nuclei. The model accurately repro- 
duces the effects observed in the experiments. Figures 17 and 18 show comparisons of 
QMSFRG predictions to experimental data for fragments from 40 Ar and 56 Fe projectiles. The 
results show good agreement between theory and experiments and indicate the larger number of 
isotopes that are produced during nuclear fragmentation. Table 3 shows the isotopic table of 
170 nuclei developed as a complete list of GCR primary and secondary nuclei to be used in 
GCR transport problems. 


20 Ne (0.6 GeV/u) + 12 C -> 2 ^ 



Zf 


Figure 5: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 20 Ne on 12 C at 0.6 GeV/u. Experimental data from Zeitlin etal. (2001). 
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We have also listed in Table 3 the mass excess (Audi and Wapstra, 1993), iso-spin, and half-life 
for the unstable nuclei along with the decay mode. This table of nuclei includes all nuclei of 
significant abundance with iso-spin +3/2 to -3 that appear in GCR transport problems. By 
contrast, early versions of the HZETRN code used a collapsed nuclear table. The expanded grid 
used here will allow for improved description of the physics, and to discuss many applications 
where unstable nuclei are central to understanding. Also listed in the first column of Table 3 is 
the index scheme that is used in the HZETRN code, which is used to label the shielding-depth 
and energy-dependent fluence matrix. The index for the nuclei are ordered by increasing mass 
number, A, followed by charge number, Z, for a given A. 


22 Ne (0.894 GeV/u)+ H ~> Z F 



Figure 6: Comparisons of the QMSFRG model to experiment for the elemental fragment 
distribution for 22 Ne on *H at 0.894 GeV/u. Experimental data from Knott eta}. (1996). 
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26 Mg (0.576 GeV/u)+ H -> Z F 



Figure 7: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 26 Mg on X H at 0.576 GeV/u. Experimental data from Knott eta/. (1996). 


27 AI (0.582 GeV/u)+ 12 C -> Z F 



Figure 8: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 27 AI on 1Z C at 0.582 GeV/u. Experimental data from Webber etaL (1990). 




28 Si (0.6 GeV/u) + 12 C -> Z F 



Figure 9: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 28 Si on 12 C at 0.6 GeV/u. Experimental data of Zeitlin etaL (2002). 


32 S + 27 AI (1 .2 GeV/u) -> Z F 



Figure 10: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 3Z S on Z7 AI at 1.2 GeV/u. Experimental data of Brechtmann etai (1988). 
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Figure 11: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 36 Ar on *H at 0.765 GeV/u. Experimental data from Knott etai. (1996). 


40 Ar (0.352 GeV/u)+ H ~> Z? 



Z F 


Figure 12: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 40 Ar on X H at 0.352 GeV/u. Experimental data from Knott etai (1996). 
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C a 0.763 GeV/u)+ H -> Z 


Expt. (Knott et al.) 
QMSFRG 


Figure 13: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 40 Ca on 12 C at 0.763 GeV/u. Experimental data from Knott etal. (1996). 


Cr 0.343 Gev/u)+ H -> Z 


Expt. (Knott et al.) 
QMSFRG 


Figure 14: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 52 Cr on J H at 0.338 GeV/u. Experimental data from Knott etal. (1996). 
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“Fe (1 .05 GeV/u)+ 12 C -> Z F 



Figure 15: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 56 Fe on 12 C at 1.05 GeV/u. Experimental data of Zeitlin etaL (1997). 


“Fe (1.05 GeV/u)+ 27 AI-> Z F 



Figure 16: Comparisons of the QMSFRG model to experiment for elemental fragment distribution 
for 56 Fe on 27 AI at 1.05 GeV/u. Experimental data of Zeitlin etal. (1997). 
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40 Ar+ 12 C (0.6 GeV/u) 



Figure 17: Comparisons of the QMSFRG model to experiment for the isotopic distribution of 
fragments for *°Ar on 12 C at 0.6 GeV/u. Experimental data from Webber etal. (1990b) . 



Figure 18: Comparisons of the QMSFRG model to experiment for the isotopic distribution of 
fragments for 56 Fe on 12 C interactions at 0.6 GeV/u. Experimental data from Webber e tal. (1990b). 
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The GCR nuclei are completely stripped and therefore the decay mode and half-life for unstable 
nuclei could be differential from those observed in laboratories on Earth. Garcia-Munoz et al. 
(1987) has noted the following differences: 1) electron capture branches, which are inactive for 
GCR nuclei relative to /J-decay, 2) for non-stripped nuclei there will be two S-shell electrons 
that participate in electron capture that will not be available in the decay of GCR nuclei, and 
3) changes in screening effects. They have made estimates of the elongation of the half-time for 
fully stripped nuclei due to these processes, which indicate an approximate doubling of the 
decay time observed for laboratory nuclei that normally decay by electron capture. Since these 
decay times are much longer than the transit time of nuclei in shielding they are not considered 
here. However, it will be useful in the future to further consider these processes when studying 
the effects of stopping GCR nuclei on planetary atmospheres or surfaces and in tissues. 



Mass Number, A 


Figure 19: Comparisons of results from the HZETRN code for the mass fluence distribution behind 
5 g/cm 2 of aluminum shielding for solar minimum conditions comparing transport with the reduced 
59-isotope grid to transport with a full 170-isotope grid. 


Figure 19 shows results from the HZETRN code at solar minimum (<j >=428 MV) behind 
5 g/cm 2 of aluminum shielding. Comparison of the mass-flux spectra for a 59-isotope grid and 
the 170-isotope grid are shown. The 170-isotope grid was developed by considering the frag- 
mentation cross-sections for a large number of GCR primary nuclei and dominant fragments in 
several materials. In Figure 19 we have scaled the fluence by the square of the ion charge as a 
measure of the ionization power of each mass group. Large differences are seen for many 
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nuclei. Figure 20 shows the percent error in the mass-fluence spectra resulting from the use of 
the reduced isotope grid for shielding depths of 5 and 20 g/cm 2 of aluminum. Errors greater 
than 100% are seen for many nuclei, however in most cases such large errors only occur for the 
less abundant nuclei. The probability of biological effects is expected to increase in a manner 
proportional to Z 2 for a given energy, and the elemental-fluence distribution may be a sufficient 
test of transport models for supporting exploration studies. In Figure 21, we show a similar 
comparison to that of Figure 20, however here for the elemental-fluence spectra. The errors are 
indeed less substantial than those of the mass-fluence spectra, yet are larger than 10% in several 
cases. Similar comparisons near solar maximum conditions (<fr=1050 MV) are shown in Fig- 
ures 22 and 23. Table 4 shows results for the elemental and neutron excess dependence of the 
point dose equivalent behind aluminum shielding. The T<-1 nuclei are not significant, while all 
other cases make important contributions to the dose equivalent. Table 5 shows the fluence at 
several depths for several of the cosmic-ray clock nuclei and other less abundant nuclei of 
interest for scientific studies. 
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Figure 20: Comparisons of the error that results from the HZETRN code for the mass fluence 
distribution near solar minimum when using a reduced 59-isotope grid compared to transport with 
a 170-isotope grid. 
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Figure 21: Comparisons of the error that results from the HZETRN code for the elemental fluence 
distribution near solar minimum when using a reduced 59-isotope grid compared to transport with 
a 170-isotope grid. 
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Mass Number, A 


Figure 22: Comparisons of the error that results from the HZETRN code for the mass fluence 
distribution near solar maximum when using a reduced 59-isotope grid compared to transport with 
a 170-isotope grid. 
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Figure 23: Comparisons of the error that results from the HZETRN code for the elemental fluence 
distribution near solar maximum when using a reduced 59-isotope grid compared to transport with 
a 170-isotope grid. 


Discussion 


I 


.t has been recognized for many years that — for the description of 
GCR transport in shielding — theoretical models and experimental data 
describing the nuclear interactions and propagation of protons, heavy 
ions, and their secondaries leading to accurate and computational efficient transport codes are 
needed. In the last 25 years, such descriptions have improved dramatically. Major milestones 
have included developing an accurate free space GCR model (Badhwar and O’Neill, 1992), the 
HZETRN code (Wilson, 1977; Wilson, et al., 1990); measuring a significant number of frag- 
mentation cross-sections (Brechtmann and Heinrech, 1988; Webber, et al., 1990b; Knott, et al. 
1996, 1997; and Zeitlin et al., 1997, 2001, 2002), and developing an accurate nuclear fragmen- 
tation model (Cucinotta, et al., 1998a). Laboratory (Schimmerling et al, 1989) and spaceflight 
(Badhwar and Cucinotta, 2000) validation data have also become available during this time 
period. The implementation of heavy ion transport models has progressed from models that 
lacked unitarity (Letaw et al., 1983) to the current fully energy-dependent models with accurate 
absorption cross-sections (Shinn et al., 1994; Wilson et al., 1993; Cucinotta, 1993). Future 
work may still be required for light-particle transport ( n , p, d, t, h, a, and mesons and their 
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decay products), including establishing production cross-section models and data, and under- 
standing the role of angular deflections, which are expected to be more important for heavy ion 
transport. However, the heavy ion problem is in much better shape, with many of the remaining 
tasks implementation issues. One exception may be improvements in fragmentation cross- 
sections and laboratory validation for the Z=1 to 5 nuclei produced from the heavier projectile 
nuclei (Z>10). 

The present paper addressed two implementation tasks: the use of a free-space GCR model, 
which includes the isotopic composition of the primaries, and the extension of the HZETRN 
code to a complete isotope grid. The problem of the isotopic distribution of the primary GCR 
and their modulation during the solar cycle has been treated in a parametric way in this paper. 

In this approach, we maintain the accuracy of the GCR modulation model for the GCR elemen- 
tal spectra from Badhwar and O’Neill (1992), however redistribute the fluence of each element 
amongst its isotopes using estimates from satellite data. Although performing more extensive 
fits to satellite data could make improvements, the coupling of the HZETRN code to the Leaky 
Box Model is suggested as the solution approach to this problem. Using the complete isotope 
grid of nuclei will allow HZETRN applications on the study of the so-called cosmic ray 
“clocks” with lifetimes similar to the time spent by GCR nuclei in the galaxy (~1 M-yr) as well 
as studies of signature nuclei from the decay of GCR nuclei in the atmosphere or on planetary 
surfaces where lifetimes on the order of a few to several thousand years will be important. 
Because the use of a reduced grid leads to error and there are no practical limitations in using a 
complete isotope grid at this time, we recommend it be used when initiatives to design space 
exploration vehicles begin. Future tasks that remain are to implement physical models of the 
GCR isotopic environment and to continue to refine the QMSFRG model, including compari- 
sons to new fragmentation data as they become available. For the many deformed and highly 
deformed projectile nuclei considered herein, methods to consider this deformation in the abra- 
sion process are needed. Another consideration for the future is the decay of radioactive iso- 
topes produced as GCR secondaries in planetary atmospheres, shielding materials, or tissue. 
Such tasks are being considered by the present authors and will be reported elsewhere. 


Acknowledgments 


T, 


he authors express their sincere appreciation for 
the partial funding support received from the projects: 
MARIE (Martian Radiation Environment Experiment on 200 1 Mars Odyssey) and NEXT 
(NASA Exploration Team), project funding numbers 111-40-10 and 732-50. 


27 



References 


Anonymous, Strategic Space Radiation Health Program Plan, NASA Office of Life Sciences and 
Microgravity Applications, Washington D.C. (1998). 

G. Audi, and A.H. Wapstra, The 1993 Atomic Mass Evaluation (I) Atomic Mass Table. Nucl. Phys. 
A565, 1-65(1993). 

G.D. Badhwar, and P.M O’Neill, An Improved Model of GCR for Space Exploration Missions. Nucl. 
Tracks Radiat. Meas. 20, 403-410 (1992). 

GD. Badhwar, and EA. Cucinotta, A Comparison on Depth Dependence of Dose and Linear Energy 
Transfer Spectra in Aluminum and Polyethylene. Radiat. Res. 153, 1-8 (2000). 

G.D. Badhwar, W. Atwell, F. F. Badavi, T. C. Yang, and T. F. Cleghom, Space Radiation Absorbed 
Dose Distribution in a Human Phantom, Radiat. Res. 157, 76-91 (2002). 

G.D. Badhwar, MArtian Radiation Environment Experiment (MARIE), Special Issue on 2001 Mars 
Odyssey, NASA-JPL Publication, in press (2003). 

C. Brechtmann, and W. Heinrich, Fragmentation Cross-sections of 32 S at 0.7, 1.2, and 200 GeV/ 
nucleon. Z.Phys. A331, 463-472 (1988). 

M. Clowdsley, J.W. Wilson, M.H.Y. Kim, R.C. Singleterry, R.K. Tripathi, J.H. Heinbockel, FT. 
Badavi, and J.L. Shinn, Neutron Environments on the Martian Surface. Phys. Medica XVII 94- 
96(2001). 

F.A. Cucinotta, GS. Khandelwal, L.W. Townsend, and J.W. Wilson, Correlations in a -a Scattering 
and Semi-Classical Optical Models. Phys. Lett. B223, 127-132 (1989). 

F.A. Cucinotta, L.W. Townsend, and J.W. Wilson, Production of 3 H at Large Momentum in a- l2 C 
Collisions at 2 A GeV. Phys. Lett. B282, 1-6 (1992a). 

F.A. Cucinotta, L.W. Townsend, and J. W. Wilson, Inclusive Inelastic Scattering of Heavy Ions in the 
Independent Particle Model. J. Phys. G: Nucl. Part. Phys. 18, 889-901 (1992b). 

F. A. Cucinotta, L. W. Townsend, and J. W. Wilson, Multiple Scattering Effects in Quasi-Elastic 
a - 4 He Scattering. Phys. Rev. C 46, 1 45 1 - 1 456 ( 1 992c). 

F.A. Cucinotta, Calculations of Cosmic-Ray Helium Transport in Shielding Materials. NASA TP- 
3354 (1993a). 

F.A. Cucinotta and R.R. Dubey, Final State Interactions and Inclusive Nuclear Collisions. NASA TP- 
3533 (1993b). 

F.A. Cucinotta, J.W. Wilson, and F.F. Badavi, Extension of BRYNTRN Code to Monoenergetic Light 
Ion Beams. NASA TP-3472 (1994a). 


28 



F.A. Cucinotta, and R. R. Dubey, Alpha Cluster Description of Excitation Energies in 12 C ( 12 C, 3a) 
X at 2. 1 GeV. Phys. Rev. C50, 979-984 (1994b). 

F.A. Cucinotta, Forward Production of Protons in Relativistic 12 C-nucleus Collisions. J. Phys. G: 
Nucl. Part. Phys. 20, 1 803- 1 8 1 5 ( 1 994c). 

F.A. Cucinotta, and J.W. Wilson, Initiation-Promotion Model of Tumor Prevalence in Mice From 
Space Radiation Exposure. Radiat. Environ. Biophys. 34, 145-149 (1995a). 

F.A. Cucinotta, Townsend, L. W., J.W. Wilson, J.L. Shinn, G.D. Badhwar, and R. R. Dubey, Light Ion 
Components of the Galactic Cosmic Rays: Nuclear Interactions and Transport. Adv. in Space 
Res. 17, 77-86 (1995b). 

F.A. Cucinotta, J.W. Wilson, M.R. Shavers, and R. Katz, The Effects of Track Structure and Cell 
Inactivation on the Calculation of Heavy Ion Mutation Rates in Mammalian Cells, Inter. J. 

Radiat. Biol. 69, 593-600 (1996a). 

F.A. Cucinotta, and J.W. Wilson, Study of Analytical Statistical Model for Decay of Light and Me- 
dium Mass Nuclei in Nuclear Fragmentation. NASA TP-3594 (1996b). 

F.A. Cucinotta, et al., Computational Procedures and Data-Base Development. In: NASA Workshop 
on Shielding Strategies for Human Space Exploration. Eds. Wilson J.W., Miller J., Konradi A., 
and Cucinotta F.A., NASA CP-3360 (1997a). 

F.A. Cucinotta, J.W. Wilson, and L.W. Townsend, Abrasion-Ablation Model for Neutron Production 
in Heavy Ion Collisions. Nucl. Phys. A. 619, 202-212 (1997b). 

F.A. Cucinotta, J.W. Wilson, J.L. Shinn, and R.K. Tripathi, Assessment and Requirements of Nuclear 
Reaction Data Bases for GCR Transport in the Atmosphere and Structures. Adv. in Space Res. 21, 
1753-1762 (1998a). 

F.A. Cucinotta, J.W. Wilson, R.K. Tripathi, and L.W. Townsend, Microscopic Fragmentation Model 
for Galactic Cosmic Ray Studies. Adv. in Space Res. 22, 533-537 (1998b). 

F.A. Cucinotta, J.W. Wilson, J.R. Williams, and J.F. Dicello, Analysis of Mir-18 Results For Physical 
and Biological Dosimetry: Radiation Shielding Effectiveness in LEO. Radiat. Meas. 132, 181- 
191 (2000). 

F.A. Cucinotta, W. Schimmerling, J.W. Wilson, L.E. Peterson, G Badhwar, P. Saganti, and J.F. 

Dicello, Space Radiation Cancer Risks And Uncertainties For Mars Missions, Radiat. Res. 156, 
682-688(2001). 

F.A. Cucinotta, Isotopic Effects in Nuclear Fragmentation and GCR Transport Problems. Proceed- 
ings of Shielding Technology Workshop, Hampton VA (2002). 

B. D. Fields, K. A. Olive, and D. N. Schramm , Cosmic Ray Models for Early Galactic Lithium, 
Beryllium and Boron Production, Fermilab-Pub-94/010-A (1994). 

M. Garcis-Munoz, J.A. Simpson, T.G Guzik, J.P. Wefel, and S.H. Margolis, Cosmic-ray Propagation 
in the Galaxy and in the Heliosphere: The Path-Length Distribution at Low Energy. Astrophys. J. 
Suppl. 64, 269-304 (1987). 


29 



A. Hesse, et al., The Isotopic Composition of Silicon and Iron in the Cosmic Radiation as Measured 
by the ALICE Experiment. 22 nd International Conference, The Dublin Institute for Advanced 
Studies, Dublin, 596-599 (1991). 

M. Kim, J.W. Wilson, R.L. Kiefer, and S.A. Thibeault, Effects of Isotope Selection on Solution 
Convergence in HZE Transport. NASA TP-3445 (1994). 

C.N. Knott, et al., Interactions of Relativistic Neon to Nickel Projectiles in Hydrogen, Elemental 
Production Cross-sections. Phys. Rev. C53, 347-357 (1996). 

C.N. Knott, et al., Interactions of Relativistic 36 Ar and 40 Ar Nuclei in Hydrogen: Isotopic Production 
Cross-sections. Physical Review C56, 398-406 (1997). 

J. Letaw, C.H. Tsao, and R. Silberberg, Matrix Methods of Cosmic Ray Propagation. Composition 

and Origin of Cosmic Rays, Maurice M. Shapiro, ed., D. Reidel Publ. Co, 337-342 (1983). 

A. Lukasiak, P. Ferrando, F.B. McDonald, and W.R. Webber, Cosmic Ray Composition of 6<Z<8 
Nuclei in the Energy Range 50-150 MeV/n by the Voyager Spacecraft during the Solar Minimum 
and Maximum Periods. 23 rd International Cosmic Ray Conference, La Jolla , USA, 539-542 
(1993). 

A. Lukasiak, F.B. McDonald, W.R. Webber, and P. Ferrando, Voyager Measurements of the Isotopic 
Composition of Sc, Ti, V, Cr, Mn, and Fe Nuclei. 24 th International Conference Ray Conference, 
Vol. 2,576-579(1995). 

K. M. Maung, John W. Norbury and David E. Kahana, Proportionally off-mass-shell equation for 

unequal mass systems, J. Phys. G: Nucl. Part. Phys. 22 No 3, 315-320 (1996). 

J.W. Norbury, F.A. Cucinotta, L.W. Townsend, and F.F. Badavi, Parameterized Cross Sections in 
Heavy-Ion Collisions. Nucl. Instrum. Meth. Phys. Res. B31, 535-537 (1988). 

E.N. Parker, The Passage of Energetic Charged Particles Through Interplanetary Space. Planet Space 
Sci. 13, 9-49 (1965). 

W. Schimmerling, J. Miller, M. Wong, M. Rapkin, J. Howard, H.G Spieler, and J. V. Blair, The 
Fragmentation of 670A MeV Neon-20 as a Function of Depth in Water. Radiat. Res. 120, 36-51 
(1989). 

J.L. Shinn, J.W. Wilson, F.F. Badavi, Fully-Energy Dependent HZETRN. NASA TP-3243 (1994). 

L. W. Townsend, J.W. Wilson, F.A. Cucinotta, and J.W. Norbury, Comparison of Optical Model 

Differences in Heavy Ion Fragmentation: Optical Versus Geometric. Phys. Rev. C34, 1491-1495 
(1986). 

R.K. Tripathi, J.W. Wilson, and F.A. Cucinotta, Medium Modified Two-Body Scattering Amplitude 
from Proton-Nucleus Total Cross-Sections. Nucl. Instrum. Meth. Phys. Res B173, 391-396 
( 2001 ). 

W.R. Webber, J.C. Kish, and D.A. Schrier, Cosmic Ray Isotope Measurements with a New Cerenkov 
Total Energy Telescope. 19 th International Cosmic Ray Conference, La Jolla , USA, 88-95 (1985). 


30 



W.R. Webber, A. Southoul, P. Ferrando, and M. Gupta, The Source Charge and Isotropic Abundances 
of Cosmic Rays with Z=9-16: A Study Using New Fragmentation Cross-Sections. The Astrophys. 
J. 348, 61 1-620 (1990a). 

W.R. Webber, J.C. Kish, and D.A. Schrier, Individual Isotopic Fragmentation Cross-Sections of 
Relativistic Nuclei in Hydrogen, Helium and Carbon Targets, Phys. Rev. C41, 547 (1990b). 

M.E. Wiedenback and D.E. Greiner, High-Resolution Observations of the Isotopic Composition of 
Carbon And Silicon in the Galactic Cosmic Rays. The Astrophys. J. 247 LI 19-L122 (1981). 

M. E. Wiedenback, The Isotopic Composition on Cosmic Ray Chlorine, 19 th International Cosmic 

Ray Conference, La Jolla , USA, 1 985, 84-87 (1985). 

J.W. Wilson, Analysis of Theory of High-Energy Transport. NASA TN D-8381 (1977). 

J.W. Wilson and F.F. Badavi, Methods of Galactic Heavy Ion Transport. Radiat. Res. 108, 231-237 
(1986). 

J.W. Wilson, L.W. Townsend, W. Schimmerling, GS. Khandelwal, F. Khan, J.E. Nealy, F.A. 

Cucinotta, L.C. Simonsen, J.L. Shinn, and J.W. Norbury, Transport Methods and Interactions for 
Space Radiations. NASA RP- 1257 (1991). 

J.W. Wilson, S.A. Thibeault, J.E. Nealy, M.Y. Kim, and R.L. Kiefer, Studies in Space Radiation 
Shield Performance. Proceedings of the Engineering & Architecture Symposium, Prairie View 
A&MUniv., 169-176(1993). 

J.W. Wilson, M.R. Shavers, F.F. Badavi, J. Miller, J.L. Shinn, and R.C. Costen, Non-perturbative 
Methods In HZE Propagation. Radiat. Res. 140, 241-248 (1994). 

J.W. Wilson, et al., NUCFRG2: An Evaluation of the Semi-empirical Nuclear Fragmentation Data- 
base. NASA TP-3533 (1995b). 

N. E. Yanasak, W.R.Binns, A.C. Cummings, E.R. Christian, J.S. George, P.L. Hink, J.Klarmann , 

R.A.Leske, M. Lijowski, R.A.Mewaldt, E.C. Stone, T.T. von Rosenvinge, and M.E. Wiedenbeck, 
Implications for Cosmic Ray Propagation from ACE Measurements of Radioactive Clock Isotope 
Abundances, Proceedings of the 26 lh International Cosmic Ray Conference, Vol. 3, 9 (1999). 

C. Zeitlin, L. Heilbronn, J. Miller, S.E. Rademacher, T. Borak, T.R. Carter, K.A. Frankel, W. 

Schimmerling, and C.E. Stronach, Heavy Fragment Production Cross-sections from 1.05 GeV/ 
nucleon 56 Fe in C, Al, Cu, Pb, and CH 2 Target. Phys. Rev. C56, 388-397 (1997). 

C. Zeitlin, A. Fukumura, L. Heilbronn, Y. Iwata, J. Miller, T. Murakami, Fragmentation Cross- 
sections of 600 MeV/nucleon 20 Ne on Elemental Targets. Phys. Rev. C64, 24902 (2001). 

C. Zeitlin, L. Heilbronn, J. Miller, A. Fukumura, Y. Iwata, and T. Murakami, Charge-changing and 
fragment production cross-sections of 28 Si on elemental targets from 400 MeV/nucleon to 1200 
MeV/nucleon, LBNL-47655 (2002). 

C. Zeitlin, et al., Results from MARIE (2003). 


31 



Table la: Comparison of HZETRN to Flight Measurements of GCR on NASA Space Shuttle and 
Russian Mir Space Station 


Mission 

DATE 

Inclination 

Altitude 

Shielding 

Dose, mGy/d 
Measured Theory 

%Difference 

Dose Eq., mSv/d 
Measured Theory 

%Difference 

STS-40 

1991 

39 

293 

Dloc2 

0.052 

0.048 

7.7 

0.13 

0.16 

-23.1 

STS-49 

1992 

28.5 

358 

Dloc2 

0.05 

0.048 

4.0 

0.127 

0.155 

-22.0 

STS-51 

1993 

28.5 

296 

Payload Bay 

0.044 

0.048 

-9.1 

0.144 

0.154 

-6.9 

STS-57 

1993 

57 

298 

Payload Bay 

0.113 

0.109 

3.5 

0.422 

0.434 

-2.8 

STS-57 

1993 

57 

298 

DLOC-2 

0.138 

0.11 

20.3 

0.414 

0.37 

10.6 

Mir-18 

1995 

51.6 

390 

P 

0.142 

0.141 

0.7 

0.461 

0.526 

-14.1 

STS-81 

1997 

51.6 

400 

0-sphere 

0.147 

0.135 

8.2 

0.479 

0.521 

-8.8 

STS-81 

1997 

51.6 

400 

Poly 3-in 

0.138 

0.138 

0.0 

0.441 

0.400 

9.3 

STS-81 

1997 

51.6 

400 

Poly 5-in 

0.129 

0.118 

8.5 

0.316 

0.368 

-16.5 

STS-81 

1997 

51.6 

400 

Poly 8-in 

0.128 

0.113 

11.7 

0.371 

0.323 

12.9 

STS-81 

1997 

51.6 

400 

Poly 12-in 

0.116 

0.111 

4.3 

0.290 

0.298 

-2.8 

STS-89 

1998 

51.6 

393 

0-sphere 

0.176 

0.148 

15.8 

0.561 

0.614 

-9.4 

STS-89 

1998 

51.6 

393 

Al 3-in 

0.167 

0.159 

4.8 

0.445 

0.488 

-9.7 

STS-89 

1998 

51.6 

393 

Al 7-in 

0.149 

0.161 

-8.1 

0.529 

0.617 

-16.6 

STS-89 

1998 

51.6 

393 

Al 9-in 

0.171 

0.162 

5.3 

0.492 

0.541 

-10.0 


Table lb: Comparisons of HZETRN Code to NASA Space Shuttle Phantom Torso Experiment on 
STS-91 (51.6 degree inclination by 390 km altitude). Measurements are taken from Badhwar etal. 
( 2002 ). 


Organ 

Measured 

(mGy) 

Theory 

(mGy) 

Theory* (mGy) 

% 

Difference 

% 

Difference* 

Brain 

2.23 

2.42 

2.26 

-8.5 

-1.4 

Bone Surface 

2.16 

2.36 

2.21 

-9.3 

-2.1 

Esophagus 

1.71 

1.79 

1.67 

-4.7 

2.2 

Lung 

1.92 

1.81 

1.69 

5.7 

11.9 

Stomach 

2.05 

2.08 

1.94 

-1.5 

5.2 

Liver 

1.88 

2.15 

2.01 

-14.4 

-6.9 

Spinal Column 

1.65 

1.98 

1.85 

-20.0 

-12.1 

Bone Marrow 

1.75 

1.98 

1.85 

-13.1 

-5.7 

Colon 

1.71 

1.9 

1.78 

-11.1 

-3.8 

Bladder 

1.58 

1.87 

1.75 

-18.4 

-10.6 

Gonad 

1.75 

1.85 

1.73 

-5.7 

1.2 

Skin/Breast 

2.46 

2.58 

2.41 

4.9 

2.0 

Skin/Abdomen 

2.35 

2.58 

2.41 

-9.8 

-2.6 


‘Comparisons that include TLD correction for LET response to heavy ions and neutrons. 
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Isotope 

6 Li* 

7 Li 

7 Be* 

9 Be 

10 Be 

10 B 

^B 

12 C 

13 C 

14 n 

15 n 

16 o 

17 o 

18 o 

20 Ne 

21 Ne 

22 Ne 

24 Mg 

25 Mg 

26 Mg 


Near-Earth Fraction 
Z^3 
0.5 
0.5 

Z=4 

0.5 

0.35 

0.15 

Z=5 

0.31 

0.69 

Z=6 

0.92 

0.08 

2=7 

0.48 

0.57 

Z=8 

0.946 

0.027 

0.027 

Z=10 

0.55 

0.10 

0.35 

Z=12 

0.64 

0.18 

0.18 


Source Fraction 

0.5 

0.5 

0.5 

0.35 

0.15 

0.2 

0.8 

0.999 

0.001 

0.78 

0.22 

0.985 

0.008 

0.007 

0.68 

0.0 

0.32 


0.74 

0.14 

0.13 






Isotope 


Source Fraction 


Near-Earth Fraction 


i i . 

2=13 


26 AI 

27 ai 


0.02 

0.98 

2=14 


0.0 

1.0 


28 Si 

29 Si 

30 Si 

32 S 

33 S 

Mg 

35 C | 

“Cl 

37 CI 


36 Ar 
l7 Ar 
l8 Ar 
°Ar 

3 Ca 

'Ca 

! Ca 

'Ca 

'Ca 


0.84 

0.08 

0.08 

2=16 

0.69 

0.15 

0.16 

2=17 

0.52 

0.41 

0.26 

2=18 

0.64 

0.03 

0.30 

0.03 

2=20 

0.4 

0.2 

0.2 

0.2 

0.2 


0.902 

0.054 

0.044 

0.96 

0.02 

0.02 

1.0 

0.0 

0.0 

1.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 

0.0 

0.0 





2=22 

0.27 

0.27 

0.31 

0.31 

0.30 

0.30 

0.09 

0.09 

0.03 

0.03 

2=23 


0.53 

0.53 

0.24 

0.24 

0.23 

0.23 

2=24 


0.26 

0.26 

0.26 

0.26 

0.48 

0.48 

2=25 


0.43 

0.43 

0.17 

0.17 

0.40 

0.40 

2=26 


0.076 

0.055 

0.084 

0.078 

0.763 

0.792 

0.076 

0.075 




Table 3 a Isotope Grid Used in Present HZETRN Calculations for Elements Z=0-8 


Index, j 

Nuclei 

z 

N 

Mass 

Excess, 

MeV 

Tz 

Lifetime 

Decay 

Mode 

1 


0 

1 

8.071323 

- 1/2 

614.6 s 

P- 

2 

'H 

1 

0 

7.28969 

+1/2 

Stable 

3 

2 H 

1 

1 

13.13572 

0 

Stable 


4 

3 h 

1 

2 

14.949794 

- 1/2 

12.33 y 

P- 

5 

3 He 

2 

1 

14.931203 

+1/2 

Stable 

6 

4 He 

2 

2 

2.424911 

0 

Stable 


7 

6 He 

2 

4 

17.5941 

-1 

0.807 s 

P- 

8 

6 Li 

3 

3 

14.0863 

0 

Stable 

9 

7 Li 

3 

4 

14.9077 

-1/2 

Stable 


11 

8 Li 

3 

5 

20.9452 

-1 

0.838 s 

P- 

13 

9 li 

3 

6 

24.9540 

-3/2 

0.178 s 

P- 

10 

7 Be 

4 

3 

15.7695 

+1/2 

53.12 d 

e 

14 

9 Be 

4 

5 

11.3477 

- 1/2 

Stable 


16 

10 Be 

4 

6 

12.6067 

-1 

1.51x10 6 y 

p- 

19 

^Be 

4 

7 

20.174 

- 3/2 

13.81 s 

p- 

12 

8 b 

5 

3 

22.92010 

+ 1 

0.770 s 

p- 

15 

9 b 

5 

4 

12.4158 

+1/2 

0.54 keV 

2ap 

17 

10 B 

5 

5 

12.0508 

0 

Stable 


20 

"B 

5 

6 

8.6680 

- 1/2 

Stable 


22 

12 b 

5 

7 

13.3689 

-1 

0.0202 s 

P- 

25 

13 b 

5 

8 

16.5623 

-3/2 

0.01736 s 

P- 

18 

io c 

6 

4 

15.6986 

+1 

19.255 s 

e+, p+ 

21 

"C 

6 

5 

10.6502 

+1/2 

0.0204 s 

e+, p+ 

23 

12 C 

6 

6 

0.0 

0 

Stable 

26 

13 c 

6 

7 

3.125011 

-1/2 

Stable 


28 

14 c 

6 

8 

3.019894 

-1 

5730 y 

P- 

31 

15 C 

6 

9 

9.8731 

-3/2 

2.449 s 

P- 

24 

12 n 

7 

5 

17.3381 

+1 

0.0111 s 

e+, p+ 

27 

13 n 

7 

6 

5.34546 

+1/2 

9.965 m 

e+, p+ 

29 

14 n 

7 

7 

2.863419 

0 

Stable 

32 

15 n 

7 

8 

0.101508 

-1/2 

Stable 


34 

16 n 

7 

9 

5.6820 

-1 

7.13 s 

P- 

36 

17 n 

7 

10 

7.781 

- 3/2 

4.173 s 

P- 

30 

14 0 

8 

6 

8.00646 

+1 

70.606 s 

e+, p+ 

33 

15 o 

8 

7 

2.8555 

+1/2 

122.24 s 

e+, p+ 

35 

16 0 

8 

8 

-4.736998 

0 

Stable 

37 

17 0 

8 

9 

-0.80900 

-1/2 

Stable 


39 

18 o 

8 

10 

-0.7821 

-1 

Stable 


42 

19 0 

8 

11 

3.3322 

-3/2 

26.91 s 

P- 

45 

20q 

8 

12 

3.7969 

-2 

13.51 s 
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Table 3h Isotope Grid Used in Present HZETRN Calculations for Elements Z=9-14 


Index, j 

Nuclei 

Z 

N 

Mass 

Excess, 

MeV 

Tz 

Lifetime 

Decay 

Mode 

38 

17 f 

9 

8 

1.95170 

+1/2 

64.49 s 

e+, p+ 

40 

18 f 

9 

9 

0.8734 

0 

0.1098 s 

e+, p+ 

43 

19 f 

9 

10 

-1.48740 

-1/2 

stable 


46 

20 f 

9 

11 

-0.01740 

-1 

11.0s 

P- 

48 

21 F 

9 

12 

-0.0476 

- 3/2 

4.158 s 

P- 

41 

18 Ne 

10 

8 

5.319 

+1 

1.672 s 

e+, p+ 

44 

19 Ne 

10 

9 

1.7510 

+1/2 

17.22 s 

e+,p+ 

47 

20 Ne 

10 

10 

-7.041929 

0 

stable 


49 

21 Ne 

10 

11 

-5.73172 

-1/2 

stable 


51 

22 Ne 

10 

12 

-8.02435 

-1 

stable 


53 

23 Ne 

10 

13 

-5.15365 

-3/2 

37.24 s 

P- 

56 

24 Ne 

10 

14 

-5.948 

-2 

3.38 m 

P- 

50 

21 Na 

11 

10 

-2.1843 

+1/2 

22.49 s 

e+, p+ 

52 

22 Na 

11 

11 

-5.1822 

0 

2.6019 y 

e+,p+ 

54 

23 Na 

11 

12 

-9.52950 

-1/2 

stable 


57 

24 Na 

11 

13 

-8.41762 

-1 

14.959 h 

P- 

59 

25 Na 

11 

14 

-9.3575 

-3/2 

59.1 s 

P- 

62 

26 Na 

11 

15 

-6.902 

-2 

1.072 s 

P- 

55 

23 Mg 

12 

11 

-5.4727 

+1/2 

11.317s 

e+, p+ 

58 

24 Mg 

12 

12 

-13.93340 

0 

stable 


60 

25 Mg 

12 

13 

-13.19275 

-1/2 

stable 


63 

26 Mg 

12 

14 

-16.21451 

-1 

stable 


66 

27 Mg 

12 

15 

-14.58654 

-3/2 

9.458 m 

P- 

69 

28 Mg 

12 

16 

-15.0188 

-2 

20.91 h 

P- 

61 

25 ai 

13 

12 

-8.9158 

+1/2 

7.183 s 

e+,p+ 

64 

26 ai 

13 

13 

-12.21032 

0 

7.17x10 s y 

e+, p+ 

67 

27 ai 

13 

14 

-17.19686 

-1/2 

stable 


70 

28 ai 

13 

15 

-16.85058 

-1 

2.241 m 

P- 

72 

29 ai 

13 

16 

-18.2155 

- 3/2 

6.56 m 

P- 

75 

30 ai 

13 

17 

-15.872 

-2 

3.60 s 

P- 

65 

26 Si 

14 

12 

-7.145 

+1 

2.234 s 

e+, P+ 

68 

27 Si 

14 

13 

-12.38503 

+1/2 

4.16 s 

e+, p+ 

71 

28 Si 

14 

14 

-21.49283 

0 

stable 


73 

29 Si 

14 

15 

-21.89506 

-1/2 

stable 


76 

30 Si 

14 

16 

-24.43292 

-1 

stable 


79 

31 Si 

14 

17 

-22.94899 

- 3/2 

0.1573 s 

P- 

82 

32 Si 

14 

18 

-24.0809 

-2 

150 y 

P- 

85 

33 Si 

14 

19 

-20.492 

- 5/2 

6.18 s 

P- 
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Table 3c: Isotope Grid Used in Present HZETRN Calculations for Elements Z= 15-19 


Index, j 

Nuclei 

z 

N 

Mass 

Excess, 

MeV 

Tz 

Lifetime 

Decay 

Mode 

74 

29p 

15 

14 

-16.9519 

+1/2 

4.140 s 

e+, P+ 

77 

30p 

15 

15 

-20.2006 

0 

2.498 m 

e+,p+ 

80 

31 p 

15 

16 

-24.44101 

-1/2 

stable 


83 

32 p 

15 

17 

-24.30534 

-1 

14.262 d 

P- 

86 

33p 

15 

18 

-26.3377 

-3/2 

25.34 d 

P- 

89 

34p 

15 

19 

-24.558 

-2 

12.43 s 

P- 

92 

35p 

15 

20 

-24.8576 

- 5/2 

47.3 s 

P- 

78 

3° s 

16 

14 

-14.063 

+1 

1.178 s 

e+, p+ 

81 

31 s 

16 

15 

-19.0449 

+1/2 

2.572 s 

e+,p+ 

84 

32 S 

16 

16 

-26.01594 

0 

Stable 


87 

33 s 

16 

17 

-26.58620 

-1/2 

Stable 


90 

M s 

16 

18 

-29.93181 

-1 

Stable 


93 

35 s 

16 

19 

-28.84633 

- 3/2 

87.32 d 

P- 

96 

Mg 

16 

20 

-30.66396 

-2 

Stable 


99 

37 S 

16 

21 

-26.89622 

- 5/2 

5.05 m 

P- 

103 

38 s 

16 

22 

-26.861 

-2 

170.3 m 

P- 

88 

33 C | 

17 

16 

-21.0035 

+1/2 

2.511 s 

e+, p+ 

91 

“Cl 

17 

17 

-24.443961 

0 

1.5264 s 

e+, p+ 

94 

35 C | 

17 

18 

-29.01351 

-1/2 

Stable 


97 

“Cl 

17 

19 

-29.52189 

-1 

3.01x10 s y 

P- 

100 

37 CI 

17 

20 

-31.76152 

- 3/2 

Stable 


104 

38 C | 

17 

21 

-29.79798 

-2 

37.34 m 

P- 

108 

39 C | 

17 

22 

-29.7998 

-5/2 

55.60 m 

P- 

95 

35 Ar 

18 

17 

-23.0482 

+1/2 

1.775 s 

e+, P+ 

98 

“Ar 

18 

18 

-30.23046 

0 

Stable 


101 

37 Ar 

18 

19 

-30.9480 

-1/2 

35.04 d 

e 

105 

38 Ar 

18 

20 

-34.7148 

-1 

Stable 


109 

39 Ar 

18 

21 

-33.242 

- 3/2 

269 y 

P- 

112 

40 Ar 

18 

22 

-35.039889 

-2 

Stable 


115 

41 Ar 

18 

23 

-33.0673 

-5/2 

109.34 m 

P- 

118 

42 Ar 

18 

24 

-34.420 

-3 

32.9 y 

P- 

102 

37 k 

19 

18 

-24.79926 

+1/2 

1.226 s 

e+, p+ 

106 

38 k 

19 

19 

-28.8017 

0 

7.636 m 

e+, p+ 

110 

39 k 

19 

20 

-33.80684 

-1/2 

Stable 


113 

40 k 

19 

21 

-33.53502 

-1 

1. 277x1 0 9 y 

P- 

116 

41 K 

19 

22 

-35.55887 

- 3/2 

Stable 


119 

42 k 

19 

23 

-35.0213 

-2 

12.36 h 

P- 

122 

43 k 

19 

24 

-36.593 

-5/2 

22.3 h 

P- 
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Table 3d; Isotope Grid Used in Present HZETRN Calculations for Elements Z=20-23 


Index, j 

Nuclei 

z 

N 

Mass 

Excess, 

MeV 

Tz 

Lifetime 

Decay 

Mode 

107 

38 Ca 

20 

18 

-22.059 

+ 1 

0.440 s 

e+, p+ 

111 

39 Ca 

20 

19 

-27.2763 

+1/2 

0.8596 s 

e+, p+ 

114 

40 Ca 

20 

20 

-34.8461 1 

0 

Stable 


117 

41 Ca 

20 

21 

-35.1375 

-1/2 

1.03x10 s y 

e 

120 

42 Ca 

20 

22 

-38.5468 

-1 

Stable 


123 

43 Ca 

20 

23 

-38.4084 

- 3/2 

Stable 


125 

“Ca 

20 

24 

-41.4691 

-2 

Stable 


128 

45 Ca 

20 

25 

-40.8125 

-5/2 

162.61 d 

p- 

131 

46 Ca 

20 

26 

-43.1350 

-3 

Stable 


121 

42 Sc 

21 

21 

-32.12109 

0 

0.681 s 

e+, p+ 

124 

43 Sc 

21 

22 

-36.1876 

- 1/2 

3.891 h 

e+, p+ 

126 

“Sc 

21 

23 

-37.8158 

-1 

3.927 h 

e+,p+ 

129 

45 Sc 

21 

24 

-41.0694 

-3/2 

stable 


132 

46 Sc 

21 

25 

-41.7587 

-2 

83.79 d 

P- 

135 

47 Sc 

21 

26 

-44.3317 

- 5/2 

3.3492 d 

P- 

138 

48 Sc 

21 

27 

-44.493 

-3 

43.67 h 

P- 

127 

“Ti 

22 

22 

-37.5483 

0 

63 y 

e 

130 

45 Tj 

22 

23 

-39.0069 

-1/2 

184.8 m 

e+, p+ 

133 

46 Ti 

22 

24 

-44.1254 

-1 

Stable 


136 

47 Ti 

22 

25 

-44.9318 

-3/2 

Stable 


139 

48 Ti 

22 

26 

-48.4871 

-2 

Stable 


142 

49 Ti 

22 

27 

-48.5581 

- 5/2 

Stable 


145 

5°Tj 

22 

28 

-51.4259 

-3 

stable 


134 

46 v 

23 

23 

-37.074 

0 

0.422 s 

e+, p+ 

137 

47 V 

23 

24 

-42.0040 

-1/2 

32.6 m 

e+, p+ 

140 

48 V 

23 

25 

-44.4747 

-1 

15.9735 d 

e+, p+ 

143 

4®V 

23 

26 

-47.9562 

- 3/2 

330 d 

e 

146 

50 V 

23 

27 

-49.2177 

-2 

1.4x10 17 y 

e+, p+ 

149 

Sly 

23 

28 

-52.1976 

- 5/2 

Stable 


152 

52 V 

23 

29 

-51.4375 

-3 

3.743 m 

P- 
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Table 3e: Isotope Grid Used in Present HZETRN Calculations for Elements Z=24-28 


Index, j 

Nuclei 

Z 

N 

Mass 

Excess, 

MeV 

Tz 

Lifetime 

Decay 

Mode 

141 

48 Cr 

24 

24 

-42.815 

0 

21.56 h 

e+, p+ 

144 

49 Cr 

24 

25 

-45.3256 

-1/2 

42.3 m 

e+, p+ 

147 

50 Cr 

24 

26 

-50.2546 

-1 

>1.8x10 17 y 

e 

150 

51 Cr 

24 

27 

-51 .4449 

-3/2 

27.7025 d 

e 

153 

52 Cr 

24 

28 

-55.4131 

-2 

Stable 


155 

s3 Cr 

24 

29 

-55.2810 

- 5/2 

Stable 


158 

^Cr 

24 

30 

-56.9287 

-3 

Stable 


148 

50 Mn 

25 

25 

-42.6216 

0 

0.284 s 

e+, p+ 

151 

51 Mn 

25 

26 

-48.2371 

- 1/2 

0.0462 s 

e+, p+ 

154 

52 Mn 

25 

27 

-50.7013 

-1 

5.591 d 

e+, p+ 

156 

53 Mn 

25 

28 

-54.6840 

-3/2 

3.74x1 0 6 

e 

159 

M Mn 

25 

29 

-55.5516 

-2 

312.3 d 

e+, P+ 

161 

55 Mn 

25 

30 

-57.7067 

- 5/2 

Stable 


157 

53 Fe 

26 

27 

-50.9414 

-1/2 

0.00851 s 

e+, p+ 

160 

M Fe 

26 

28 

-56.2485 

-1 

Stable 


162 

55 Fe 

26 

29 

- 57.4751 

-3/2 

2.73 y 

e 

164 

“Fe 

26 

30 

-60.6013 

-2 

Stable 


167 

57 Fe 

26 

31 

-60.1760 

-5/2 

Stable 


163 

55 Co 

27 

28 

-54.0239 

-1/2 

17.53 h 

e+,p+ 

165 

“Co 

27 

29 

-56.0353 

-1 

77.27 d 

e+, p+ 

168 

57 Co 

27 

30 

-59.3400 

- 3/2 

271.79 d 

e 

166 

“Ni 

28 

28 

-53.900 

0 

6.077 d 

e+, p+ 

169 

57 Ni 

28 

29 

-56.0757 

-1/2 

35.60 h 

e+, p+ 

170 

»Ni 

28 

30 

-60.2233 

-1 

Stable 
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Table 4a: Elemental (Z) and Neutron Excess (Y) Dependence on GCR Dose Equivalent (cSV/yr) 
Behind 0 g/cm 2 of Aluminum Shielding 








Table 4b: Elemental (Z) and Neutron Excess (Y) Dependence on GCR Dose Equivalent (cSv/yr) 
Behind 5 g/cm 2 of Aluminum Shielding 







GCR Pose Equivalent per Year near Solar Minimum 

Z Y<0 Y=0 Y=1 Y=2 Y=3 Y>3 Total-Z 











Table 5: Annual Fluence of Several Unstable GCR Nuclei Versus Depth of Aluminum Shielding 
Near Solar Minimum (^>=428 MV) 



Fluence per cm 2 per year 

Nuclei 

0, g/cntf 

2, g/crrt 

5, g/cirt 

20, g/cirt 

6 He 

0.0 

9-lxlO 1 

2.1x102 

5.6x102 

8 Li 

0.0 

I.OxlO 3 

2.2x103 

5.4x103 

9 Li 

0.0 

3.0x1 0 1 

7.2x10i 

1.8x102 

10 Be 

8.8x103 

9.1x103 

9.3x103 

8.8x103 

14C 

0.0 

7.1x102 

1.7x103 

3.4x103 

180 

9.9x103 

9.3x103 

8.4x103 

5.1x103 

26AI 

5.7x1 0 2 

7.6x102 

9.4x102 

1.0x103 

36CI 

8.1X10 2 

7.7x102 

7.1x103 

4.4x102 
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